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Abstract 

■ We consider a quantization of the Bianchi IX cosmological model based on 



taking the constraint to be a self-adjoint operator in an auxiliary Hilbert 
space. Using a WKB-style self-consistent approximation, the constraint cho- 



> 

o 

Q\ • sen is shown to have only continuous spectrum at zero. Nevertheless, the 

o. 

auxiliary space induces an inner product on the zero-eigenvalue generalized 
eigenstates such that the resulting physical Hilbert space has countably infi- 

O" 1 

nite dimension. In addition, a complete set of gauge-invariant operators on 

the physical space is constructed by integrating differential forms over the 

spacetime. The behavior of these operators indicates that this quantization 

preserves Wald's classical result that the Bianchi IX spacetimes expand to a 

maximum volume and then recollapse. 
ICP:0460 
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I. INTRODUCTION 



Finite dimensional cosmological models are a mainstay of modern general relativity for 
the simple reason that they address many exciting conceptual issues without forcing the 
researcher to confront the complicated, infinite dimensional, and in the quantum context 
perturbatively non-renormalizable formalism of full 3+1 gravity. Our interest here will be 
the quantum theory and our motivation for studying this model is exactly that stated above. 
Again, we ask the historical questions (e.g., |ffl-ffl): "What are the quantum observables?" 
"How are they constructed or even defined?" And, "On what Hilbert space, if any, do they 
act?" 

A proposal for addressing these questions was introduced in [|J and applied to finite 
dimensional parameterized Newtonian systems, the free relativistic particle, and "separable, 
semi-bound" models (which include the Kantowski-Sachs and locally rotationally symmetric 
(LRS) Bianchi IX cosmological models, see [§J). For such systems, the approach was found 
to define a Hilbert space of physical states, to produce a complete set of gauge-invariant 
quantum operators, and to reproduce (in a certain sense) the classical feature that the semi- 
bound models describe spacetimes with homogeneous spatial hypersurfaces first expand and 
then recollapse. The strategy was based on that of DeWitt || and takes seriously the idea 
that observables in General Relativity are given by integrals of differential forms over the 
spacetime manifold. The analogues of such integrals were, in fact, performed in the quantum 
theory. 

In order to implement this basic proposal, a mathematical framework must first be 
introduced. The main framework of || is based on that of Dirac [[| and considers an 
auxiliary Hilbert space in which gauge-dependent operators act. For the models in question, 
a physical inner product was induced on the generalized eigenvectors of the constraint by 
this auxiliary structure despite the fact that no solution of the constraint was normalizable 
in the auxiliary inner product. 

However, j4| considered only relatively simple models. While some models addressed by 
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section V of || are not explicitly solvable even classically, the usual minisuperspace models 
(Bianchi I, Kantowski-Sachs, and LRS Bianchi IX) that fall into this category are solvable 
and have even been deparameterized |7[. Also, the separable nature of these systems was 
explicitly used in || to construct the physical Hilbert space and to study convergence of the 
integrals that define the gauge invariant operators. 

Our goal here is to apply these same ideas to the more complicated Bianchi IX, or 
Mixmaster, minisuperspace l^-jlOfl. Like other Bianchi models, it describes a class of solutions 
to the vacuum Einstein equations that are foliated by a family of homogeneous spacelike 
hypersurfaces. This preferred family effectively reduces the gauge symmetry of full General 
Relativity to that of time reparameterizations. 

Unlike the simpler models, the Bianchi IX cosmology has not been completely solved or 



deparametrized classically. In fact, it is often |TT],|T2j considered to be an example of chaotic 
behavior in General Relativity. It therefore constitutes a much more stringent test of the 
proposal of ||. We shall see that, by introducing new tools and a WKB-style self-consistent 
approximation, all of the corresponding results are obtained. 

We begin our presentation by reviewing the method of || in section [TT]. Section pTI| 
then applies these techniques to the mixmaster model, beginning with a description of the 
quantum theory (section |III A| ) and followed by a study of solutions to our constraint QUI B|) 
and a characterization of the physical Hilbert space ( [III Q ). One consequence of our results 
is that, within the validity to this approximation, we show that the space of solutions to 
the Hamiltonian constraint for Bianchi IX is infinite dimensional. Finally, a complete set 
of gauge-invariant operators is constructed and the recollapsing behavior of our quantum 
theory is verified in [HI D| . We close with a summary of the results and a comparison with 
other work. 
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II. THE GENERAL APPROACH 



Because our goal is to apply the techniques of j| to the Bianchi IX minisuperspace, 
we now summarize the relevant points of B for the reader's convenience. The general 
setting considered is that of a time reparameterization invariant system with Hamiltonian 
constraint: 

h = (2.1) 

for some function h on the phase space T = R 27V (see [ffl for a more general treatment). 
Here, as in ||, we will take capital letters to stand for quantum operators and lower case 
latin letters to refer to their classical counterparts. For each coordinate qi,Pj on T, a 
family Qi(t), Pj(t) of self-adjoint operators is introduced on an (auxiliary) Hilbert space 

n = l 2 (r n ). 

Now, choose some time independent symmetric factor ordering H{Q i {t) ) Pj(t)) corre- 
sponding to h(qi,pj) and consider the associated self-adjoint operator H. The families 
Qi(t), Pj(t) are chosen to satisfy 

HA 

— = i [N(t)H,A] (2.2) 

for any A = A(Qi(t), Pj(t)) and a family of lapse operators N(t) that are proportional to 
the identity on 7i; N(t) = n(t)ll. The commutators of the families form a quantized version 



of a generalized Peierls bracket [13[] and the formulation of [§]] provides a notion of gauge 
transformations corresponding to time reparameterizations. 

This structure is to be applied toward two goals: i) The construction of a physical 
Hilbert space Tiphys and ii) The construction of a complete set of gauge-invariant operators 
that commute with H . Here, we use the same notion of a complete set as in |4j]. That is, 
we say that a set of operators on TL p h ys is complete if the corresponding set of phase space 
functions obtained in the classical limit separates the points of the phase space. For param- 
eterized Newtonian systems and so-called "separable semi-bound" systems (which include 
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the Kantowski-Sachs and LRS Bianchi IX cosmological models), both of these objectives are 
attained in J3|. We will now see that they may be accomplished for the full-fledged Bianchi 
IX model as well. 

The physical Hilbert space is to be defined as follows. We first find a convenient "energy 
basis" of states \n,E) for E G [A, 0] (for some A < 0) and n G Z + that satisfy H\n,E) = 
E\n, E) and 

(n',E'\n,E) = 5(E-E')5 rhn , (2.3) 

and span the subspace of TC corresponding to the interval [A, 0] in the continuous spectrum 
of H. The physical space TC p h ys is to be the closure of the vector space generated by the 
symbols {|n, 0)} for n G Z + with the physical inner product: (\n', 0), \n, 0)) p hys = 8n,n'- 

For a gauge-invariant operator A on TC that commutes with H, we would then like to 
consider an operator A p h ys on TC p h ys given by: 

A phys \n,0) = J2°<nk{0)\k,0) (2.4) 
k 

where the coefficients a n k are defined by A\n, E) = J2k a nk(E)\k, E). This procedure is 
satisfactory when a n k{E) is continuous in E on [A, 0] and in this case A p h ys is a symmetric 
bilinear form on TC p h ys when A is symmetric on TC. In this way our physical inner product 



captures the classical reality conditions |14| to the same extent as the auxiliary inner product 
on TC. 

Such gauge invariant operators were constructed in ||] by integrating over the parameter 
t that labels the families of operators that appear in |2.2| . In particular, given two operators 
A and Z on TC, we follow and consider 

/CO f) 
JtA(t)(-e(Z(t)-r)) (2.5) 

and 

/OO f) 
Jt{-9{Z{t)-r))A{t) (2.6) 
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where 0(Z(t) — r) is the projection onto the non-negative spectrum of Z(t) — r. The object 

[A]z=T = \{[A}z=r + [A] R z=T ) (2.7) 

is formally symmetric and in the classical limit is simply a function on the space S of solutions 
that, when applied to a particular solution s, is a linear combination of the values of A at 



the points along s at which Z — r. Furthermore, when the integrals ^5] and ^6] converge, 
they define operators that commute with H For appropriate systems and choices of A 
and Z, this procedure was shown in || to produce a complete set of well-defined operators 
([A}z= T ) P hys on Hphys- Because the correspondence between the gauge invariant A and the 
physical operator A p h ys is direct, we now drop the subscript phys from these operators. 

One final property derived in || is that a quantization of this type captures the classical 
"recollapsing behavior" of separable semi-bound cosmological models. That is, in contrast 
to quantization schemes based on the Klein-Gordon inner product there is a quantum 
analogue of the classical result that these models expand and then recollapse. For such 
models, if Z is taken to be the scale factor a = ln(det g) for the metric g on a homogeneous 
spatial hypersurface, then the operators [A] a=T vanish in the large r limit. Because the 
corresponding classical quantities satisfy [a] a=T — > as a direct consequence of the fact that 
a given solution s will not reach arbitrarily large values of a, we interpret 

[A] a=T -> as r -> +oo (2.8) 

as the quantum statement of recollapse. The Bianchi IX model also recollapses classically 



16fl and we shall see that, when quantized as above, it recollapses quantum mechanically as 



well. 



III. THE QUANTIZATION OF BIANCHI IX 

The discussion of [f|] considered the parameterized Newtonian particle and separable 
semi-bound models, but our interest here is in the more complicated Bianchi IX minisu- 
perspace. We will see that the same (qualitative) results can be obtained as we will now 

6 



construct i) a non-trivial physical Hilbert space from generalized eigenvectors of the con- 
straint operator as well as ii) a complete set of gauge invariant operators and Hi) verify that 
our quantum theory leads to the recollapsing behavior derived classically in [[RJ . 

We follow much the same path as in section V of 0]. Section A| presents the system 



and sections [Tl fB] and |111 C| construct a convenient "energy" basis for the auxiliary space 7i 



in terms of which the physical states can be simply characterized. This part of the analysis 
is carried out at the somewhat physical level of WKB-style self-consistent approximations. 
Section D fD] considers operators of the form [A] Q=T for a = ln(detg) and g the three- 



metric on such a slice and finds that matrix elements of these operators converge for a 
complete set of appropriate A so such [A] a=T exist as bilinear forms on Ti p h ys - We then modify 
the construction of [A] Q=T by introducing a "level operator" L whose physical interpretation 
is less clear, but which improves the behavior of our operators. The modified [^4] a=r can 
then be proven to define (for appropriate A) a complete set of bounded operators on TL p h ys - 
Both the forms [A] Q=T and the operators [A] Q=T vanish as r —>■ oo, demonstrating that our 
model preserves the classical recollapsing behavior. 

A. The Quantum System 

The mixmaster model describes anisotropic homogeneous solutions to the 3 + 1 vacuum 



Einstein equations. For this case, the metric g^- on a homogeneous slice |10| may be written 
in the form 

9H = e 2a (e 2 % (3.1) 

for a matrix given by diag{(3 + + Vo/3~, (3 + — \/3P~, — 2(3 + ) which defines the anisotropy 
parameters /3 and the scale factor a. Taking the metric to be in this form fixes all of the 
gauge freedom except that of time reparameterizations. 

Since the classical phase space is described by (a, (3~ ,p a ,Pp+,P/3~) £ R 6 we take the 
auxiliary space TL to be L 2 (R 3 ,c/ 3 x) following Q. Our quantum theory will be built from 
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the configuration operators a,^ which act on functions if) E L 2 (R 3 ) by multiplication as 
well as the derivative operators P a = —ih-J^, Pp± = —ih-r^ which are the corresponding 



momenta. 



Recall that the full Hamiltonian constraint for this model may be written 10 



as 



H = -Pi + P| + + P|_ + e* a (V((3) - 1) (3.2) 

using a rescaled lapse N(t) = (^ E ) 1//2 e~ 3a iV(t) and Hamiltonian H = (^r) 1//2 e 3a ff. (where 
N and H are the usual ADM lapse and Hamiltonian of general relativity jl7| applied to the 



minisuperspace ansatz p.l|) . Note that by writing |3T2j as a quantum operator, we have chosen 
a certain factor ordering of the classical expression. The explicit form of the potential is 
given by 

V(J3) = \e~^ + - \e~ w+ cosh 2y/3f3~ + 1 + ^e 4/3+ (cosh 4^/3" - 1) (3.3) 

which has a global minimum V — at the origin f3 ± — but has no other critical points. 
The operator [T2]is symmetric on L 2 (R 3 ) and we shall assume that some suitable extension 
is in fact self-adjoint. This extension will also be written as H and no distinction will be 
made between the two operators. 



The operator H of |0| is Hermitian and we use it to introduce a "time dependent" set 



of configuration bases {\a, (3 ; t')} which are related through 

\a,f3 + ,f3-;t') = e i ^'- £ > H \a,f3 + ,f3-;t). (3.4) 

The (time dependent) direct product decomposition H = H a]t <8> H^t induced by the fac- 
torization \a, P + , P~; t) = \a;t) <g> \/3 + ,/3~;t) = \a; t) |/? + , f3~; t) will be used heavily in what 
follows. 



As our system is time reparametrization invariant and |3^ is classically constrained to 
vanish, it should be emphasized that this t merely serves as a label and is not the Newtonian 
time parameter of non-relativistic quantum mechanics. Nevertheless, it behaves similarly 
when used as a technical tool and we use it to construct the families Qi(t) and Pj{t). These 
will be defined by Q;(0) = Q h Pj{0) = P 3 and 
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Qi(t) = exp (iH J dt'N{t'])Qi[ -iH J dt'N(t') 
Pj(t) = exp (iH J* dt'N(t'j)Pi( - iH J* dt'N(t'j) (3.5) 

for some N(t) = n(t)lL, n(t) G R. This "lapse operator" and its relation to gauge transfor- 
mations are discussed in [HI. 



The generalized eigenvectors of the constraint |3.2| with eigenvalue zero will form the 
physical Hilbert space TC p hy S - In order to study solutions of H\ip) = 0, we will consider the 
operator 

H X = P| + + P|- + e^{V{(3) - 1) (3.6) 

and for each t', introduce a family H*(a; t') of operators on ?ip-t> parametrized by a: 

H x = [ da\a;t')(a;t'\® H*(a;t'). (3.7) 

This H*(a;t') is a Hermitian operator on the Hilbert space TCftt' = L 2 (R 2 ) of square inte- 
grable functions of f3 + and f3~ . We similarly define P^± and /3 . Note that, for various a, 
the H*(a) differ only in the overall scale of the potential V r *(/3 ± *; a) = e 4a (l/(/3 ± *) — 1) and 
that, because the region in which the potential is less than E for any £eR has finite area 
(see appendix A), the spectrum of these operators is entirely discrete |[8| and the degeneracy 
of each spectral projection is finite. 

Thus, we may label the eigenstates of H*(a;t') in Tip-f as \n;a;t')/3 for n G Z + . The 
state \n;a;t')p is to have the eigenvalue E n (a) (which does not depend on if) of H*(a;t') 
and these states are orthonormal: p(n'; a; t'\n; a; t')p = 5„, n '- Since the family H*(a;t') is 



analytic in a, the states |n; a; t) may also be chosen smooth in the label a for each n [19 
and we may express any state mTC as 



J da^2f n (a)\a)\n;a)p (3.8) 



for some f n (a). Until D ID], we will work at some fixed value of t' and suppress this label on 



the states, operators, and spaces. 
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B. Solutions of the Constraint 



In the following, we operate at the somewhat physical level of self-consistent approxi- 



mations and to use 3.8 to study the eigenvectors of H; that is, to solve H\ip) = E\ip). We 



begin by considering the action of H on 3.S: 

#|V> = J daJ2(h 2 Q^(fn(a)\n;a)p) + f n (a)E n (a)\n;a) )\a) (3.9) 
Now, from the equation 

{H*{a)-E n {a))\n;a) p = (3.10) 
in Tip it follows that the states \n; a)/3 may be chosen to satisfy 

h H] ^ = (H*(a)-E n (a))J - Ae4a{V ~ ^ ^ 

= M mn \m; a) p (3.11) 

\m\ a 



, "7/3 



where ^ H *^} E (7m]_ ls ^ ne inverse of H*(a) — E n (a) except when acting on states 
with E m (a) = E n (a) (i.e., states degenerate with \n;a)/3), which it annihilates. 

The matrix M mn is antihermitian since M nm vanishes when E n = E m and, for E n ^ E m , 
we have 

M mn = (m;a\ ~ 4 e 4a (V - l)|n; a) 

(H*(a) - E n (a))± 

= (m;a\-4e^(V-l)\n; a) 
E m (a) - E n (a) 

for which the numerator is Hermitian and the denominator is real and antisymmetric. We 



now use M„ m to rewrite the derivative in 3.S as 



d d 

-fr(52fn{u)\n\a)p) = + ^2fm{a)M nm )\n;a)i 3 . (3.13) 

(Jul n n uu. m 

Our goal is to study generalized eigenvectors of H for which the solutions \ip) of [3]8] are 
at least delta- function normalizable. Thus, for almost every a the list of functions f n (a) 
should give the components (/(«))„ of an a-dependent vector f(a) in the Hilbert space 
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I 2 of square summable sequences that satisfies J da\f(a)\ 2 = (ip\ift). On this space, the 
coefficients M mn define an antihermitian operator through (Mf) n = Y^ m M nrn f ' m . Let U(a) 
be the corresponding unitary path ordered exponential Vexp(— f£ da'M(a')) satisfying 
"' - U~ X M and let / = U~g for some g. Then we have 



da 



—f n + Y. M nmf m ) = (U^g) n (3.14) 



and we may write |3T^ as 

H\i>) = J daJ2{h 2 (U—g) n + (Ug) n E n (a))\n; a). (3.15) 

It follows that \ip) is an eigenvector of H with eigenvalue E exactly when the corresponding 
g is a generalized L 2 function and satisfies 

2 

h 2 —g+(E(a)og) = Eg. (3.16) 

Here, E(a) = U~ 1 (a)E(a)U(a) and E(a) is the operator on I 2 such that (E(a)f) n = 
E n (a)f n . We will show that for E < 0, there is one such solution for every g 6 I 2 . 

The equation |3.16| to be solved for g(a) is essentially a one-dimensional time independent 



Schrodinger equation for a multiple component wavefunction with potential — E(a) and 
energy — E. Thus, existence, normalizability, and multiplicity of the solutions to |3.16| is 
related to whether this "potential" confines a (fictitious) particle to finite values of a or lets 
such a particle escape to the asymptotic regions. Note that, as the potential V*(j3*;a) of 
H*(a) becomes more and more shallow for large negative a, E n (a) —>■ as a — > — oo and our 
fictitious particle is unbound in this direction for E < 0. In particular, for large negative a, 
even the ground state energy of H*(a) is positive so that —E n {a) approaches from below 



as a — > — oo and there can be no normalizable eigenvector with E < 0; all solutions to |3.16 
with E < are delta-function normalizable in the region a < 0. 



Solutions of |3.16| with E > are difficult to analyze in the a — > — oo limit without a 
uniform bound on the rate at which E n (a) vanishes. Note, however, that the existence of 
such a unifofm bound would show that the spectrum of H is discrete for E > so that zero 
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is the largest continuum eigenvalue in dramatic contrast to the separable cases of Q. A 
WKB approximation of the type used below for a — > +00 might provide insight but is even 
more difficult to control than at the positive a end. As we only require knowledge of the 
spectrum of H on some [a, b] containing zero, we take a < and 6 = and do not concern 
ourselves with the case E > 0. 

The remaining task is to study the behavior of E n (a) for large, positive a. To do so, we 
write g in the form 

g = Vexp[f a da'N(a')]g Q (3.17) 
J 

where V = Pexp[/ a da'N(a')] is the path ordered exponential that satisfies |^ = NV. Our 



equation |3.16| can now be written as 

[N 2 V + NV + EV]g = EVg (3.18) 

where N = dN/da. In order to show that a solution exists for all g , we will (approximately) 
solve the matrix equation 

N 2 V + NV + EV — EV = (3.19) 



As in the finite dimensional WKB method, assume that iV 2 >> N, so that 3.19 becomes 
iV 2 = EH. — E. Since E is a positive self-adjoint operator with known spectrum, a number 



of square-roots y EH. — E are easily defined. For the moment, let iV be any one of these. 



Now, recall that E(a) = U^ 1 (a)E{a)U(a). It follows that N = U- 1 (a){^ EH - E)U(a) 
and we have 

9fTM < dU rr-llrrr,. tt™ 



— (UV)-[—U- 1 }UV = U— = ^EH.-EUV (3.20) 
oa da da 



which has the solution V = U~ X V exp[/ a (y EH — E)—M\. This form of V will be convenient 
for the approximations below. 

We now reason as follows. Consider again H*(a) and the associated eigenvalue problem 

E n ip n = -h 2 Mn + e Aa {V - 1)W (3-21) 
12 



Here, A is the two-dimensional Laplacian and ip n is the wavefunction ip n = (/3 + , [3~\n). Since 



3.21 may also be written as 



n h 



2 ^ n = -A^„ + -y(V-l)Vn (3.22) 



the ratio E n /h 2 must depend only on % and a through e 4a /h 2 . Similarly, we must have 
e~ 4a E n = f(e 4a /h 2 ) for some function / and it follows that a — > oo is essentially just the 
model's classical limitQ However, since the potential (V — l)e 4a has only a single critical 
point (which is, in fact, a global minimum) we expect that as h — > every energy level falls 
toward this minimum. Thus, e~ Aa E n — > — 1 as a — > oo. 

This will in turn imply that, for large a, the nth state is concentrated near (3 — 0. Since 
E n = (n\H*\n), —A > 0, and (V — 1) > —1 as a — ► oo, the behavior of E n (a) implies that 
(n| (1/ — 1) \n) — > —1. But (V — 1) = —1 only at the origin, so we must have J- R 2_ Bc \ip n \ 2 
for every e-ball B t centered at ^ = 0. 

For small enough e, the potential V — 1 within B e is essentially that of a (rotationally 
symmetric) Harmonic oscillator (V ~ 8e 4a [/? +2 + /?~ 2 ]). Thus, since E n — > — e 4a for every n, 
for each > there is some a for which if* (a) may be replaced by the rotationally invariant 
Harmonic oscillator Hamiltonian —A + 8e 4a (/3 2 _ + Z? 2 .) — e 4a when considering states with 
n < N to some accuracy e . At this point, it is convenient to replace the label n with the two 
occupation numbers n\,n% of the oscillator so that E nitU2 (a) ~ —e 4:a + 2^/2e 2a h(ni + n2 + l). 
Similarly, we find (m 1 ,m 2 \(V- l)\n x ,n 2 ) w c {mum2 ) {nun2 )e- 2a , for c (mi)m2 ) (ni)n2 ) independent 
of a and vanishing for |mi — n\\ > 2 or |m2 — n.2 1 > 2, so that 



M ^ C L (mi,m 2 )(ni,n 2 ) oq\ 

JW (mi,m 2 )(ni,n 2 ) ~ 2rvfe^ i \ yo.Z,0 



which is also independent of a in this regime. On the other hand, y EH — E — > ue 2a for 
some u 2 = II so that we may neglect M in comparison to this term and write 

V w C/ _1 exp(ue 2a ). (3.24) 



Thanks to Charles W. Misner for pointing this out. 
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Note in particular that N 2 ~ e 4a while N ~ e 2a so that for large a we do indeed have 
N 2 » N. 

C. Eigenvectors and the Physical Hilbert Space 



While |3.24| (approximately) solves |3.19| , the question arises of which solutions V(a)g are 



(delta-function) normalizable. Suppose first that u has some eigenvector ~g with eigenvalue 
+ 1. Then, the corresponding solution g = Wg has diverging I 2 norm (|#| 2 ) as a — > oo. 
However, for any eigenvector g Q of u with eigenvalue —1, the corresponding g = Wg Q is 
square integrable on the region a > (J a>Q \g\ 2 < oo) for all square summable g Q . We thus 
take u — — 1 for the general normalizable solution. 

For large negative a, we already know that the behavior of our solutions is that of a free 



particle. It therefore follows that there is one delta-function normalizable solution to |3.16 
for each g e I 2 and every E < so that the continuous spectrum of H includes the interval 
(— oo, 0]. This verifies that H is of the type described in section |T| and provides one of the 
labels (E) of our useful basis of states. 

In fact, we have shown that all states in 7i are of the form 

mUE)}) = [ dE\g (E),E) = [ dE^2(g (E)) n \Xn,E) (3.25) 

where \g Q ,E) is the solution of H\ip) = E\ip) given by g(a) = VE(a)g Q for g e I 2 and 
(Xn)m = 8 n ,m- If, for convenience, we now refer to the states \x n , E) as \n, E) then we have 
found a basis for 7i such that H\n,E) = E\n,E) and (n, E\n' , E') = 6(E — E')5 n>n /. The 
physical Hilbert space TC p h ys may now be constructed as in Q as the closure of the span of 
{\n, 0)} with respect to the inner product (|n, 0), \m, 0)) p h ys = <5n,m- 

D. Physical Operators 



Having constructed Ti. p hy S} we now address both the formation of physical operators on 
this space as in 0] through integration over t and the verification of recollapsing behavior. 
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We first construct a complete set of operators that i) has a straightforward physical interpre- 
tation in the classical limit, ii) are defined as bilinear forms on the physical space, and Hi) 
satisfy the recollapse criterion |2.8| as bilinear forms. We will then regularize these operators 
using a new operator L whose physical interpretation is less than clear. The regularized 
operators, however, form a complete set of bounded operators on the physical Hilbert space 
which satisfy in the sense of strong convergence; that is, convergence on a dense set of 
states |2( 

Consider then the objects [A\ Q=T of ^L7\ . From J|, the action of these operators on \n, 0) 
may be written as 

(|n*, 0), [A] L a=T \n, 0)) phys = 2m(n*, 0\A[H, 9(a - r)} |n, 0) (3.26) 

and 

(|n*, 0), [A]* =T \n, 0)) phys = 2m(n\ 0| [H, 9(a - r))A\n, 0) (3.27) 

Now, for any bounded operator B on TC, let 

B re9 = tA-#(« " r')B6(a - r')-^—. (3.28) 
ri + % id — i 

This B reg is bounded since H is self-adjoint and is Hermitian when B is. As in the 
operator [B reg ] a=T has the same classical limit as [-B] a=T for r' < r since, when the constraint 
is satisfied, the factors of combine to have no effect and the extra step functions 9(a — r') 
are classically irrelevant for r > r'. 

To see that the [B reg ] a=T give a complete set of operators, note that B may be chosen to 
be any bounded function of (3 or of Pp± . By also inserting into B appropriate functions of 
l±sign(P a ) as in [|J, we may construct objects [B reg ] a=T whose classical limits give arbitrary 
(bounded) functions of /3 ± or Pp± evaluated at point where a = t and either P a > or 



P a < 0. However, from the result of JT6[ each classical solution has a = r only twice (if at 
all), once with P a > and once with P a < 0. Thus, since we consider the [B reg ] a=T for all 
real r, this is a complete set of operators for the Bianchi IX model. 
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Because all of the solutions to |3.16| that were used to build TL p h ys are square-integrable 



over the region a > (see section |111 C| ), the states 6(a — r')\n*,0) and 8(a — r)\n,0) are 
normalizable in 7i. It follows that the physical matrix elements (|n*,0), [B reg ] a=T \n, 0)) p h ys 
are well-defined and finite since, for example, 

(|n*,0>, [B reg ] L a=T \n } Q)) phys = 27r(n*,O|0(a - r')\B9(a - T ')-^—]9{a - r)\n,0) (3.29) 

and B9[a — t)j^- is bounded on 7i. This means that integrals defining [B reg ] a=T converge as 
bilinear forms on TC p hy S - Also, since \8(a — r)\n, 0)| 2 — > as r — > oo, these objects converge 
to zero as forms and verify the classical recollapsing behavior. 

Now, a modification of |2.5| and |2l^ leads to objects that can be proven to be bounded 
operators on TC p hy S - To determine if some [A] a=T is a well-defined operator, it must be 
checked whether the norms of the states obtained by acting with [A] Q=T on physical states 
are finite. Thus, from |3.26| and |3.27| , we are interested in sums of the form 



^2\(n*,0;t'\O\n,0;t')\ (3.30) 



n* 



for O = [A, H]8(a — r) and O = 8(a — r)[if, A]. Equivalently, we can replace [A, H] in |3.30 
with AH or HA as appropriate since H vanishes on the physical subspace. 

In order to improve the convergence of |3.30| , we consider the new "level operator" L = 
L © 11 ± where the direct sum corresponds to the decomposition H = H © 7i L of TL into Ti 
and its orthogonal complement. The operator IL^ is the identity operator on Ti, 1 - and L is 
defined by 

L \n, E) = n\n, E) (3.31) 



so that L is self-adjoint. Now, instead of 2.5 and 2.6, consider the integrals 



[A)Lr = r f(L(t))A(t)p(a(t)-r)f(L(t)) 



X' 



— R r°° a 

[A] a=T = J ^ f(L(t))-9(a(t) - r)A(t)f(L(t)) (3.32) 

where / is an as yet unspecified function that vanishes for large L. Here, L(t) = L since 
[H, L] = from |3.31| , so that |3.32| amounts to multiplying expressions |0| and |2l] by f(L) 
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~L R 

on the left and on the right. As before, L4] q=t ' = [A] Q=r for self-adjoint A, [A] a=T = 

|(P3L + P3L)^nd 

(\n*, 0), [A] L a=r \n, 0)) phys = -i2it(n\ 0\f(L)AH9(a - r)f(L)\n, 0) 

(In*, 0), PfL-ln, 0}) phys = i27r(n\ 0\f(L)6(a - r)HAf(L)\n, 0) (3.33) 



Again, let B reg = jf^d(oi — t')B9(ol — t')-^^. For this case, the matrix elements in |3.33 
may be written 

(\n*,0),[BZ g ] L a=T \n,0)) phys = -2nf(n*)f(n)^ n .(r')\Be(a - t')-^—\Mt)) (3.34) 

where |^„(t)) = 0(a(t') — r)\n,0) is a state in Tt with finite norm z T {n). Let us choose 
/ such that \f(n)\ < nz and t' such that r' < r. Thus, / may depend on r' but not 
on r. Since z T (n) is a nonincreasing function of r, we have ,z T (n) < z T i(n), and the matrix 
elements |3.34| are bounded in absolute value by 



\\B\\ z (n) 

2vr \f(n)f(n*)\ ^(71)^(^)11511 < iLJl^lon (3.35) 

nn* z T i[n) 

where \ \B\ \ is the operator norm of B on 7i. 

The physical norms of [A]^.|n, 0) are then bounded by 

R It (n*) 

ll[^e 9 ]_|n,0)||^ s < E^^(2vr) 2 (3.36) 

and [B reg ] a=T is a bounded operator on H.ph ys . Furthermore, the action of [B reg ] a=T on any 



of the states \n, 0) vanishes as r — ► oo since both bounds in |3.36| become zero in this limit. 
This is clear for the first bound and follows for the second from the facts that the sum 
over n* converges and that z T (n*) is a positive nonincreasing function of r that vanishes as 
t — » oo. It follows that [B reg ] a converges strongly to zero and, in fact, [B reg ] a=r converges 
uniformly to zero. 

Again, in the classical limit and on the constraint surface, the factors of A cancel and 
the factors of 6(a — r') have no effect for r' < r. In addition, since [L, H] = 0, L is conserved 
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along a classical solution and, for a given solution s, [f(l)j^9(a — T')b9(a — T')j^f(l)] a=T (s) 
tends to zero for large r if and only if [b] a=T — > as r — > oo. We thus conclude that our 
model describes a quantum theory of recollapsing cosmologies. 

Since I is some function on the phase space, it can be expressed in terms of /3 ± and pp± . 
But b may be an arbitrary bounded function of (3^ and pp± so that {[/(O] 2 ^} is still sufficient 
to separate points on the phase space. In this way, we have constructed a complete set of 
quantum operators for the Bianchi IX minisuperspace. 



IV. DISCUSSION 

As claimed in the introduction, we have applied the techniques of [|J to the Bianchi IX 
cosmological model and thereby construct a quantum theory. Specifically, a Hilbert space 
structure on the generalized eigenvalues of H was induced by the auxiliary space TC aU x, 
a complete set of observables was constructed, and recollapsing behavior was derived, all 
despite the complicated and perhaps chaotic nature of the model. 

Finally, since the literature on the mixmaster cosmology is extensive, a short comparison 
of section |T| with previous work is presented below. Note first that the observables studied 
here are, in spirit, much like those of pT| (and suggested by Moncrief), which describe 
anisotropies and momenta at the maximal value of a. Such objects are easily written down 



in the formalism used here by choosing Z in 2.7 and 2.6 to be P a and r to be zero, since P c 



vanishes classically on the maximal volume surface. Note that |21] also defines observables 
away from maximal volume which, for volumes suitably close to maximal are (exactly) 
unitarily related. This "unitary evolution" then breaks down at a finite separation from 
maximal volume. In our case, because the physical operators [A] a=T vanish for large r, we 
also find that the [A] Q=T and [A] a=r / are not unitarily related when the separation between 
r and r' is large. On the other hand, there is no reason to think that any of the [A] a=T of 
PTI| are unitarily related or that there should be a preferred finite separation at which such 
a unitary evolution breaks down. Thus, we expect these quantizations to be inequivalent. 
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We now single out the recent construction of exact solutions to the quantum constraint 
by Moncrief and Ryan |22[] for more detailed discussion. Such solutions were found earlier 



by Kodama for the constraint written in terms of Ashtekar variables, but the work of 
f22|| uses the variables of section JTJ so that the comparison with our approach is direct. A 
similar comment applies to the super-symmetric solutions of Graham 

However, the quantum states of |22[ do not solve the constraint p.2| as they are related 



to a different factor ordering of the classical expression. Instead, they solve 

o = *V = |r* - b§-j - ^ - ^ + +><y - i» (4.i) 



for some real B where, following [22], we have dropped the Ks. 

Thus, our setting does not coincide with that of |22] and some effort will be needed to 
reconcile them. In particular, |22] does not consider an auxiliary Hilbert space and we must 



introduce one here for comparison with section |T|. Note that [^2| provides explicit solutions 



only for B = —6 7^ and that for B e R, B 7^ 0, the differential operator above is not 
Hermitian in the Hilbert space L 2 (R 3 , dadf3 + df3~) . This operator is, however, self-adjoint in 
the Hilbert space L 2 (R 3 , e~ Ba dadj3 + d[3~), which we will therefore take to be our auxiliary 
space. 

A rescaling of the wavef unction ip = e Ba ^ 2 (p leads to the equivalent constraint 

- *** - 55* " " Wf* ~ (JV + el " [V ~ * (4 ' 2) 



on G L 2 (R 3 ,<i 3 x). Note that the arguments of |Hj apply as well to [4.2| as to |3.2j , the 
only difference being a relative shift of E n (a) by Thus, the behavior at large a is 
unchanged and all energies are "bound" at this end. The difference appears at large negative 
a, where the behavior of our wavefunction reduces to that of a free particle. If H"<p = E(j), 
then this free particle has energy —(E+ ^-). Thus, the (probable) upper bound on the 
continuous spectrum moves from zero to — ^j-. In particular, the solutions described by 



22| vanish rapidly as a — > +00 but become asymptotically constant as a — » —00. The 
rescaled wavefunction <p then behaves as e~ Ba ^ 2 for large negative a so that such solutions 
are normalizable whenever B < such as when B = — 6. 
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Thus, zero is now in the discrete spectrum of the constraint and the construction of pTC 



is unnecessary since the relevant eigenvectors form a subspace of the auxiliary space and 
so inherit the inner product directly. Also, because the spectrum is discrete, from 0] we 
expect that the integrals |2.5| and [2.6| which attempt to define our gauge-invariant operators 
will now fail to converge. 

Finally, if we allow B to be imaginary then the original operator H' is self-adjoint on 
L 2 (R 2 ,d 3 x). The same analysis follows as for B real, but now the boundary between the 
discrete and continuous spectrum moves to E = \B\ 2 /4 > so that zero again lies in the 
continuous spectrum. Solutions of the form described by |22| would then correspond to 



elements of H p h ys of |Tl| and the construction of gauge-invariant operators proceeds as in 
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APPENDIX A: FINITENESS OF THE REGION V < E 

We now quickly derive the fact that the area A(E) of the region of the /3 + ,/3~ plane 
on which V of |3l| is less than E is finite for any E in R. To do so, we use the 27r/3 
discrete rotational symmetry of V and its well-known "triangular" shape Recall that 
this "triangle" points along the (5 + axis. 

For large positive let h E (j3 + ) be the value of /3_ for which V{(3 +) (3J) = E. In this 
part of the plane, there are always two solutions ±h E (j3 + ) and we take h E {(3 + ) to be the 
positive one. A(E) is then finite if we have 
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oo 

E 



A 



h h ((3 + )d(3 + < oo 



for any finite A. 

From 13.31, it follows that 



V(P+,P-) < -\e~ w + + ^e 4/3 +(cosh4v / 3/?_ - 1) 
so that h E ((3 + ) is less than the h E (f3 + ) defined by 

= -2e~ 2/3+ + e 4/3+ (cosh4v / 3^( / 9+) - 1) 
Clearly, h E — > for large (3 + so that, if (3 + is large enough, cosh4\/3/i 



It follows that h E <h E < |e~ 6/3 + and that [£]] holds. 
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